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SUMMARY 


Existing frequency domain modal frequency and damping analysis methods 
are discussed. The effects of truncation in the Laplace and Fourier trans^ 
form data analysis methods are described in detail. Methods for eliminating 
truncation errors from measured damping are presented. Implications of 
truncation effects in fast Fourier transform analysis are discussed. Limited 
comparison with test data Is presented. 
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INTRODUCTION 


Flight flutter testing is generally a time-consuming procedure. It 
involves the installation of complex excitation generators such as vanes, 
inertia exciters, or impulsive devices (ref, l) as well as the response 
transducers and the associated electronic equipment. During flight testing, 
many flights are required to fully explore the aircraft critical flight spec- 
trum, producing a large amount of test data which must be' subsequently analyzed. 

Considerable simplification in equipment Installation may be obtained if 
turbulence excitation can be used instead of mechanical excitation. In any 
event, atmospheric turbulence and buffet degrade the response data from all 
types of mechanical excitation, except for random excitation, where it would 
most probably help more than hinder (ref. 2). Thus, the availability of 
suitable random response analysis methods, in addition to the existing har- 
monic analysis methods, would be a great advantage. The random analysis 
methods, like the current harmonic analysis methods, place the burden of data 
reduction on the computer, which, when used in the interactive mode with the 
test engineer, can provide a basis for real-time flutter testing. 

The exciter installation and data acquisition and analysis problems are 
further compounded in space shuttle type vehicles , where weight is of para- 
mount importance and the cost of exploring the entire critical flight spectrum 
with many flights prohibitive. The nonstationary nature of the flight envi- 
ronment and the relatively short duration of each flight within the atmosphere 
place a premium on the need for transmitting as much response data as possible. 
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and as quickly as possible ^ to the groiand station. An increase either in the 
rate of sampling of the transducers , or in the number of transducers ^ is pos- 
sible, if the data sample length can be reduced without a loss in the accuracy 
of the analysis . 

A reduction in the sample length of the random response data is accom-- 
panied by a reduction in the statistical accuracy of the frequency domain 
modal response spectra. The statistical accuracy can be restored at the 
expense of resolution through a corresponding Increase in the effective anal- 
ysis bandwidth. This increase in analysis bandwidth produces a truncation 
effect in the response spectra. The truncation effect can occur in the fre- 
quency domain modal analysis derived from the Fourier transform of not only 
the impulse response time history but also the cross- and auto-correlation 
functions of response due to random and impulse-type excitations. 

The effect of truncation is studied by using a single-degree-of-freedom 
system. Existing frequency domain harmonic analysis methods are briefly 
discussed as an introduction to the truncation effect and to illustrate the 
format of the data presentation. 


HARMONIC ANALYSIS 


The simplest method for obtaining aircraft modal frequency and damping 
data is through stick pulse generated free decay data (figs. 1 and 2). How- 
ever, narrow band filtering is required both to isolate each mode in turn and 
to minimize noise due to the presence of turbulence. Computerized least 
squares fit methods such as the Moving Block Analysis (ref. 3) can be used 
to obtain damping data from the log decrement of the decay once the resonant 
frequencies have been identified by spectral analysis. 

Stick pulses, in general, may not excite all the modes of interest and 
may produce an uneonservative estimate of the damping. For close resonances, 
narrow band filtering may not isolate each mode, resulting in a beating decay 
response (fig. l). Under such circumstances, it is possible to extract 
meaningful data only if the modal damping and amplitudes are comparable in 
each of the modes. It is, however, possible, through the Fourier transform, 
to transform the decay data into the frequency domain (figs, 3 and and 
thereby resolve the modes. 

This Fourier transform process can be illustrated mathematically by 
considering the relationship between the response y(t) of a linear system and 
a general force x(t), given by 


y(t) 



h(r ) X (t-r) dr 


( 1 ) 
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where h(T) is the impulse response function of the system at time T. For a 
single-degree-of-freedom system, the impulse response function is given by 

-6o) T I ' " ' 

h(r) = — ' Sinw /l-6^ (2) 

/ 2 ^ 
m o) 1 r- 6 
r V 

where : 

m is the generalized mass, 

co^ is the ang\ilar resonant frequency, and 

5 is the viscous damping coefficient. 

If the force is of sufficiently short duration that it can be considered 
to be an impulse I6(t), where 5(t) is the delta function, then the response 
time history reduces to 


y(t) = h(t) 1 (3) 

The response spectrum y(ia)) , obtained from Fourier transform of the 
time history (eq. (l)), is related to the force spectrum x(iw) by 

y(ico) = H(ioj) x(io)) (U) 

where H(iw) is the frequency response function of the system. For a single- 
degree-of-freedom system, 


m 


/ 2 2 
(w - CO 

r 


2i6co^co). 


(5) 


The Fourier transforms of the response y(t) and the force x(t) are 
defined by 


and 


y(iw) 


x(ico) 


ij£y(t) e-^.t 


( 6 ) 


(T) 
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respectively. Again^ if x(t) can be considered an impulse 16(1)9 then the 
force spectrxm reduces to 


x { io)) 


I 

27T 


( 8 ) 


and the response spectrum to 


yCicj) 


H(ig;) I 
27T 


( 9 ) 


■which is simply the system frequency response function multiplied by a 
constant . 

Two formats can be used in the presentation of the frequency domain 
response data. In the response modulus vs frequency presentation (fig. 3)9 
the resonant frequency is located approximately at the peak response 9 and the 
viscous damping coefficient 9 which is twice the structure damping, is obtained 
by dividing the half power point bandwidth by twice the resonant frequency. 

It is accurate for well-separated modes . For close modes, as the modiolus 
represents the total response vector from the origin and not necessarily the 
modal vector, errors in the measured modal frequency and the viscous damping 
coefficient may result. The extraction of modal damping may even be prevented 
by the failure to resolve the half power points (fig. 3). 

To overcome these limitations, both amplitude and phase are retained and 
presented in a format of a Nyquist plot (fig. U) in which the real part of 
the response is in phase and the imaginary part is oub of phase relative to 
the excitation. This method of modal analysis was first suggested by Kennedy 
and Pancu (ref. U) . The resonant frequency is located at the point on the 
curve where the rate of change in arc length with frequency is at a maximum. 
The viscous damping coefficient Is obtained from the half power points as 
previously described or by first measuring the angle 0 subtended at the modal 
origin, by the arc between any frequency point f and the resonant frequency 
point f^9 and then using the relationship 

f - f 

<5 = ^ cot(0) (10) 

A strong feat-ure of the Nyquist plot response data representation is 
that mode shapes can be identified by means of the modal response vector. 

The more common method of generating the response Nyquist plots is by 
means of a slow sine sweep using mechanical in-flight excitation, such as 
inertia exciters, in which the force output is used as reference. This method 
has been computerized for multimodal analysis (ref. 5) ? employing a least 
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squares curve fit technique to minimize the effect of extraneous noise and 
used in a computer/test engineer interactive mode for flutter testing. 

In transforming the free decay data as previously discussed, no trunca- 
tion effects -were ohserved due to the relatively high modal damping and the 
need for including one heat as a minimimi in the decay sample. In the second 
example of a stick pulse excited decay (fig, 2), the decay was prematurely 
truncated after one and five seconds to illustrate the effect on the frequency 
domain response (fig. 5)* The Nyquist plots of the response become more oval 
in appearance as the decay sample duration is progressively reduced. If the 
Nyquist plots are analyzed by the conventional method described above, uncon- 
servative estimates of the damping are obtained, (See table 1.) In order 
to obtain useful damping data from these Nyquist plots, a method eliminating 
the effect of truncation from the damping must first be developed. 


TRUNCATION THEORY 


Due to the similarity between the cross-correlation and the impulse- 
response functions with the auto-correlation function Ryy ( r) of a single- 
degree-of-freedom system excited by a constant spectrum force, Sp (ref, 2) 
and defined by 


Ryy(r) 


TrSp e 
_ — 

2m cj ^6 
r 


[Cos oj 



r 



Sin CO 

r 



( 11 ) 


it is only necessary to describe the equations for any one of the above 
functions. The impulse-response function and the cross-correlation function 
of a single-degree-of-freedom system, when excited by constant spectrum force, 
exist only for positive time. 

If the Laplace transform or the single-sided Fourier transform of the 
autocorrelation function of the response Ryy (r) is used, with limits of inte- 
gration from zero to infinity, instead of the full Fourier transform, phase 
information is retained in the response spectrum (ref. 2). The resulting 
response spectrum S(ia)) is given by 

S / . V 

S(iw) = - H(ico) + 25 H(ia))) (12) 

Umo) 6 \^r / 

r 


The characteristic response function 


S(io)) 


has properties similar to the 


frequency response function E(ia)). This method provides a powerful tool in 
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modal response analysis of random as well as impulse response data in the 
frequency domain. The previously described methods for extracting the modal 
damping and frequency can be employed as long as no truncation effect is 
present . 

If we let y(t) represent any of the above time functions, with the 
understanding that they exist only for positive time, the complex response 
spectrum y(ico) is obtained from Fourier transform of y(t). 


y(ico) 


^ oo 

' y(t) 

. oo 

-iojt 
e at 

where 

y(t) = 

0 

for t < 0 


( 13 ) 


In reality, the response time history is truncated at some finite 
time Tm. The resulting estimated response spectrum y(ia)j^) (ref . 6) is 
given by the relationship 


A. V X 

y(ic^3_) = 2^ 


/ 

•/o 


f 

2tr J 

•'—ot 


ni -ico t 

y(t) e ^ at 


” -ico t 

D(t) y(t) e at 


(lU) 


where D{t) is the weighting or the truncation function. 

Three weighting functions (ref. T) are consiaerea in this paper. They 
are the "ao-nothing" or the boxcar weighting, generally aefinea by 


D(t) 


1 for-Tui<t<Tm 

0 elsewhere 


( 15 ) 


the Hanning weighting function defined by 

D(t) = -| (l + Cos for - T^<t<r^ 

= 0 elsewhere 

and the Bartlett weighting function defined by 

D(t) = (l - for-T <t<Tj^ 

= 0 elsewhere 


(l6) 


(IT) 
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The reverse Fourier transform for the response Is given by 


y(t) 



oo 

j(iw) e 


iort; 


dw 


(18) 


On substituting for y(t) in equation (lU) and rearranging the order of the 
integration, the estimate of the spectrum becomes 


where 


y(iw) 



-i(o; -w)t 

e dt y(iw)dw 



OO 

Q (w^-w) y(ico) 


dw 


Q(co^-cu) 



-i(o; -w)t 
I dt 


( 19 ) 


( 20 ) 


Q(o) -co) is referred to as the spectral -window. The weighting functions 

defined by equations (l6) to (l8) and the corresponding spectral windows are 

illustrated in figure 6. In the application discussed in this paper, the 

spectral windows are complex (ref, 6 and 8) since the weighting functions are 

one sided and exist only in positive time from zero to r . 

m 


TRUNCATED DATA REDUCTION 


For a linear system excited by random force (or impiilse) of constant 
spectral density, the response spectrum y(io;) is proportional to the frequency 
response function of the system. Equation ( 19)5 with y(ico) replaced by the 
frequency response function of a single-degree-of-freedom system and a 
constant force spectrimi, has been integrated by using contour integration 
for the ”do~nothing” and the Bartlett weighting in references .6 and 8, respec- 
tively. It has recently been solved by the author for the Hanning weighting. 

A typical effect of the truncation due to the Hanning weighting is illustrated 
in figure T* The single-degree-of-freedom response plots have been normalized 
relative to the untruncated plot. The other two weighting functions differ 
only in the degree of truncation effect. The ^^do-nothing” weighting function, 
while exhibiting the smallest truncation effect, suffers from the imdesir able 
side lobes (fig. 5) which may be mistaken for modes or may interfere with 
other nearby modes. The Bartlett weighting function suffers a greater resolu- 
tion loss, as can be seen by comparing figure 8 with figure 5- 

The resonant frequency is still identified by the peak rate of change of 
arc length with frequency, but the procedure for estimating the damping from 
the truncated curves is different from the methods previously described. At 
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first (ref. 6 and 8)5 the damping coefficient was extracted with the assist- 
ance of a nondimens ional parameter defined hy the peak rate of change of arc 
length with frequency^ divided by the radius of curvature at the resonant 
frequency 5 from theoretically predicted curves. In these curves 5 the above 
parameter is plotted as a function of the resonant frequency multiplied by 
the true damping coefficient. These curves were originally developed for 
use in high frequency structural response studies and consequently are 
unsuitable for flutter data analysis due to the relatively low aircraft 
response frequencies. 

A more useful graphical format 5 which provides direct damping readout , 
is presented in this paper and illustrated in figure 9 for the Hanning trun- 
cation. The measured damping coefficient 6* is plotted against the true 
damping coefficient 6 as a function of the ratio of the effective data 
analysis bandwidth Af divided by the resonant frequency. The effective 
analysis bandwidth Af is related to the maximimi truncation time by 


Af = 



* 

The measured damping coefficient 6 is defined by 


(21) 



(22) 


where p is the radius of curvature at resonance, and 

ds is the arc length at resonance contained within 
the frequency interval of df 

It can be observed that as the maximum truncation time becomes large, 
the measured viscous damping coefficient approaches the true value . 

This method of obtaining the damping from the truncation-affected single- 
degree-of-freedom system Hyquist plots has been computerized, for potential 
use in real-time analysis. The number of iterations required to converge to 
the correct damping from the estimated damping is illustrated in figure 10. 

The convergence is carried out in two or three sequences and is very rapid. 

The number of steps in the initial sequence is selected to speed up the 
iteration, especially in cases of severe truncation. 

The damping of the free decay record (fig. 2 ) as obtained by the comput- 
erized method for the "do-nothing”, Hanning, and Bartlett truncations, a 
least squares fit to the free decay, and the restored Nyquist plot method 
(fig. 11) are summarized in table 2 . The method of restoring the Nyquist 
(ref. 9) plot involves the weighting of the decay with a known exponential 
weighting to meet the required 55 dB dynamic range criteria (ref, lO) for 
the decay. Analysis is thereafter carried out conventionally and the damping 
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corresponding to the exponential weighting is subtracted from the measured 
damping to arrive at the true modal damping. It is more common to apply the 
exponential weighting function to the correlation function. This method has 
been used in flight flutter testing in England (ref. 11). The results from 
the analysis of the one’-second decay record indicate the existence of a 
lower bound on the accuracy for the above frequency domain analysis methods. 


TRUNCATION IN POWER SPECTRAL DENSITY 


A method based on the cross-spectral analysis previously discussed was 
developed to predict the truncation effect in power spectral density (PSD) 
analysis. The effect of the truncation on the normalized PSD is illustrated 
in figure 12 for the Hanning weighting. A computer program was developed to 
obtain the damping from the 3 dB points by using the quadratic .curve fit. 

A graphical method for obtaining the true damping coefficient 6 from the 

A 

measured damping coefficient 6 for various ratios of effective analysis 
bandwidth to resonant frequency is illustrated in figure 13. 

A Hanning smoothed power spectral density plot of a typical aircraft 
response to high-speed buffet is illustrated in figure ik , Due to the very 
high speed, no reliable flutter test data are available for comparison. The 
analysis bandwidth of 0,5 Hz produces a truncation error in the two predom- 
inant modes at ID. 2 Hz and lU.6 Hz, On allowing for the truncation effect, 
the viscous damping coefficients from the measured 3 dB point values of 0.11 
and 0.0^^ are reduced to 0.068 and 0.02 for the two frequencies, respectively . 
This method suffers from the same inaccuracies as the modulus method pre- 
viously discussed. It does, however, provide an indication of the damping 
where none previously existed. 


FAST FOURIER TRANSFORM AND TRUNCATION 


The above methods have been basically developed for the Blackman and 
Tuckey type of analysis (ref. 7)* Truncation effects occur also in the fast 
Fourier transform (FFT) method of analysis. An indication as to the type of 
truncation present in FFT analysis of cross spectra is obtained from refer- 
ence 12. The expected cross-spectral estimate E [Sxy(f,T,k)] is given by 

E rSxy(f,T,k)] = Rxy(r) (23) 

As the cross-correlation function of a single-degree-of-freedom system 
excited by white noise is one sided, as previously discussed, it is concluded 
that the estimated cross spectrum obtained from FFT analysis is subject to 
Bartlett truncation errors. 
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The effect of truncation on the normalized PSD and cross spectral peak 
response is illustrated in figure 15 as a function of the resonant frequency 
multiplied by the maximum time delay and the viscous damping coefficient, 
fr m 5 observed that for the ’^do-nothing^* truncation, the curve 

reaches unity near 1* This corresponds not only to the damping 

criteria for cross-spectral analysis established in reference 9? but also to 
the rule of thumb for PSD resolution for the analysis bandwidth to be one- 
fourth of the 3 dB point response bandwidth. 

Attention is drawn to the fact that the Bartlett truncation curve con- 
verges to unity very slowly. Thus the use of cross-correlation functions 
obtained from the indirect method of first computing the spectra using the 
FFT and then transforming to time domain, will not only have the Bartlett 
truncation error but also an additional truncation error in transforming 
from the frequency domain to the time domain. These truncation errors in 
correlation functions are discussed in references 13 » 1^ and 15. Thus a very 
large number of transformation points must be used in determining the correla- 
tion function through the indirect method. 


CONCLUSION 


Methods for eliminating truncation errors from modal frequency and 
damping data have been presented for the cross- and power-spectral analysis. 
These methods have the potential for use in a computer/test engineer inter- 
active mode, for random and impulsive-type response data analysis. Future 
work will include an evaluation of the methods against simulated and real 
flutter test data with buffet and turbulence excitation and the study of 
truncation effects in FFT-type analysis involving multiple Fourier transform 
operations • 
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TABLE 1. VISCOUS DAMPING COEFFICIENT OF TRUNCATED 
PLOTS DETERMINED BY USING EQUATION (lO). 


Decay 
Time , 

"Do-Nothing" 

Hanning 

Bartlett 

Seconds 

Truncation 

Truncation 

Truncation 

1 

0 . l86 

0.336 

0.2h8 

5 

O.OhQ 

0.0T3 

0.068 


TABLE 2. COMPARISON OF VISCOUS DAMPING 
COEFFICIENT BY VARIOUS METHODS. 


Decay- 
Time 5 
Seconds 

Least 

Square 

Decay 

Restored 

Nyquist 

Plot 

TRUNCATION THEORY 

"Do-Nothing" 

Hanning 

Bartlett 

1 

O.Oi+5 

i 

0.092 

0.037 

- 

0.035- 

0.059 

5 ^ 

^ ! 

0 . 038 

1 

O.OhO 

0.030 

0.037 

O.OI+I+ 
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Figure 1. Typical Filtered Stick Ptilse Excited Delay with Beating. 



Stick Pulse Excited Decay. 





Figure k, Fourier Transform of Beating Stick Pulse Excited Decay 





FREQUENCY Hz 



Fig-ure 5- Fourier Transform of Stick Pvilse Excited Single 
Mode Decay Showing Truncation Effect.. 





Figure 6. Spectral Windows and Weighting Functions. 
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NORMALIZED g 

CROSS SPECTRUM 

df 



= 0.015 


= 0.003 


IMAGINARY 
PART OF 
NORMALIZED 
CROSS 
SPECTRUM 


Figure T. The Effect of Truncation on the Normalized Cross 
Spectrum of a Single-Begreei-of-Freedom System Excited 
hy White Noise Hanning Weighting* 


FREQUENCY Hz 



Figure 8, Fourier Transform of Stick Pulse Excited Decay with 

Bartlett Weighting. 



ACTUAL DAMPING - S 


Figure 9« Method for Obtaining Damping from Trimcation 
Affected Cross Spectra of a Single-Degree-of-Freedom 
System - Hanning Weighting. 

































by High-Speed Buffet , 





NORMALIZED PEAK PSD AND CROSS SPECTRUM 



Figure 15 . Effect of Truncation on the Peak PSD and Cross Spectral 
Response Resolution of a Single-Degree-of-Freedom System 
Excited by White Noise, 
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